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Topology and geometry are essential to our understanding of modern physics, underlying many foundational 
concepts from high energy theories, quantum information, and condensed matter physics. In condensed matter 
systems, a wide range of phenomena stem from the geometry of the band eigenstates, which is encoded in 
the matrix-valued Wilson line for general multi-band systems. Using an ultracold gas of Rb atoms loaded in 
a honeycomb optical lattice, we realize strong-force dynamics in Bloch bands that are described by Wilson 
lines and observe an evolution in the band populations that directly reveals the band geometry. Our technique 
enables a full determination of band eigenstates. Berry curvature, and topological invariants, including single- 
and multi-band Chern and Z 2 numbers. 


Geometric concepts play an increasingly important role in 
elucidating the behavior of condensed matter systems. In 
band structures without degeneracies, the geometric phase 
acquired by a quantum state during adiabatic evolution el¬ 
egantly describes a spectrum of phenomena [1]. This geo¬ 
metric phase—known as the Berry phase—is used to formu¬ 
late the Chern number [2], which is the topological invari¬ 
ant characterizing the integer quantum Hall effect [3]. How¬ 
ever, condensed matter properties that are determined by mul¬ 
tiple bands with degeneracies, such as in topological insula¬ 
tors [4, 5] and graphene [6], can seldom be understood with 
standard Berry phases. Recent work has shown that such sys¬ 
tems can instead be described using Wilson lines [7-10]. 

Wilson lines encode the geometry of degenerate states [11], 
providing indispensable information for the ongoing effort to 
identify the topological structure of bands. For example, the 
eigenvalues of Wilson-Zak loops (i.e., Wilson lines closed 
by a reciprocal lattice vector) can be used to formulate the 
Z 2 invariant of topological insulators [7] and identify topo¬ 
logical orders protected by lattice symmetries [8, 9]. Al¬ 
though experiments have accessed the geometry of isolated 
bands through various methods including transport measure¬ 
ments [3, 12, 13], interferometry [14, 15], and angle-resolved 
photoemission spectroscopy [5, 16], Wilson lines have thus 
far remained a theoretical construct [7-10]. 

Using ultracold atoms in a graphene-like honeycomb lat¬ 
tice, we demonstrate that Wilson lines can be accessed and 
used as versatile probes of band structure geometry. Whereas 
the Berry phase merely multiplies a state by a phase factor, 
the Wilson line is a matrix-valued operator that can mix state 
populations [11] (Fig. lA). We measure the Wilson line by 
detecting changes in the band populations [17] under the in¬ 
fluence of an external force, which transports atoms through 
reciprocal space [18]. In the presence of a force F, atoms 
with initial quasimomentum q(0) evolve to quasimomentum 



FIG. 1. Wilson lines and effectively degenerate Bloch bands. 

(A) In a non-degenerate system (left), adiabatic evolution of a state 
through parameter space R results in the acquisition of a geomet¬ 
ric phase factor, known as the Berry phase. In a degenerate system 
(right), the evolution is instead governed by a matrix-valued quantity 
called the Wilson line. If the degenerate levels can be experimentally 
distinguished (blue and yellow colouring), then population changes 
between the levels are detectable. (B) The band structure of the low¬ 
est two bands of the honeycomb lattice in effective energy units of 
|F| d, where F is the applied force used to transport the atoms and 
d is the distance between nearest-neighbor lattice sites. As |F| is 
increased, the largest energy scale of the bands becomes small com¬ 
pared to |F| d. At large forces (iii), the effect of the band energies 
is negligible and the system is effectively degenerate. In this regime, 
the evolution is governed by the Wilson line operator. We distin¬ 
guish between the bands using a band mapping technique that detects 
changes in the band population along the Wilson line path. 


q(t) = q(0) + Ff/fi after a time t. If the force is sufficiently 
weak and the bands are non-degenerate, the system will un- 
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dergo adiabatic Bloch oscillations and remain in the lowest 
band [18]. In this case, the quantum state merely acquires a 
phase factor comprised of the geometric Berry phase and a 
dynamical phase. At stronger forces, however, transitions to 
other bands occur and the state evolves into a superposition 
over several bands. 

When the force is infinite with respect to a chosen set of 
bands, the effect of the dispersion vanishes and the bands 
appear as effectively degenerate (Fig. IB). The system then 
evolves according to the formalism of Wilczek and Zee for 
adiabatic motion in a degenerate system [11]. The unitary 
time-evolution operator describing the dynamics is the Wil¬ 
son line matrix [19]: 

Wq( 0 )-,q(t) = Pexp[i dqAq], (1) 

where the path-ordered (V) integral runs over the path C in 
reciprocal space from q(0) to q(t) and Aq is the Wilczek- 
Zee connection, which encodes the local geometric properties 
of the state space. 

In a lattice system with Bloch states |T>q) = 
in the n* band at quasimomentum q, where f is the 
position operator, the elements of the Wilczek-Zee con¬ 
nection are determined by the cell-periodic part \u^) as 

Aq’’^ _ i(rt^|Vq|rtq ). The diagonal elements (n = n') are 
the Berry connections of the individual Bloch bands, which 
yield the Berry phase when integrated along a closed path. 
The off-diagonal elements (n 7 ^ n') are the inter-band Berry 
connections, which couple the bands and induce inter-band 
transitions. 

Although the evolution described by Eq. 1 must be path- 
ordered when the Wilczek-Zee connections at different quasi¬ 
momenta do not commute, it can also be path-independent 
under certain circumstances [9, 20]. For example, when the 
relevant bands span the same Hilbert space for all quasimo¬ 
menta, as is the case in our system, the Wilson line operator 
describing transport of a Bloch state from Q to q reduces to 
Wq^q = [9, 19, 21]. Consequently, the elements 

of the Wilson line operator simply measure the overlap be¬ 
tween the cell-periodic Bloch functions at the initial and final 
quasimomenta [9, 21]: 

= «|«^). ( 2 ) 

Hence, access to the Wilson line elements facilitates the char¬ 
acterization of band structure topology in both path-dependent 
and path-independent evolution. In both cases, the topological 
information is encoded in the eigenvalues of the Wilson-Zak 
loops. In the latter case, the simplified form of the Wilson line 
in Eq. 2 additionally enables a map of the cell-periodic Bloch 
functions over the entire BZ in the basis of the states \uq) at 
the initial quasimomentum Q. 

We create the honeycomb lattice by interfering three blue- 
detuned laser beams at 120(1)° angles (Fig. 2A). At a lattice 
depth Vb=5.2(l)F^^, where F^^=h^/( 2 mA|) is the recoil en¬ 
ergy, Xl is the laser wavelength, and m is the mass of ^^Rb, 
the combined width e ^ h x 3 kHz of the lowest two bands 
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FIG. 2. Reaching the Wilson line regime in the honeycomb lattice. 

(A) Schematic of the honeycomb lattice in real space with A (B) 
sublattice sites denoted by solid (open) circles. A static lattice is 
formed by interfering three in-plane laser beams (blue arrows) with 
frequency u. Sweeping the frequency of beam i by Sui creates a 
force Fi in the lattice frame in the propagation direction of beam 
i [19]. (B) Two copies of the first BZ of the honeycomb lattice, 
separated by a reciprocal lattice vector G. By changing the relative 
strengths of Fi (red arrows), the atoms can be moved along arbitrary 
paths in reciprocal space. Each BZ features non-equivalent Dirac 
points K and K' at the corners of the hexagonal cell. High-symmetry 
points r, at the center of the BZ, and M, at the edge of the BZ, 
are also shown. (C) The population remaining in the first band for 
different forces after transport to T -h 0.2G (green), T -h 0.55G (red), 
and r -h G (blue). Inset numbers i to iii refer to band schematics 
in Fig. IB, representing the diminishing effect of the dispersion for 
increasing force. The data agree well with a two-level, tight-binding 
model (dashed line) which approaches the Wilson line regime (thick 
shaded line) at large forces. Discrepancies at larger forces result from 
transfer to higher bands and match well with ab initio theory using a 
full band structure calculation including the first six bands (thin solid 
line) [19]. For all subsequent data, we use \F\d/e = 4.8, indicated 
by the dashed gray line. Error bars indicate the standard error of the 
mean from ten shots per data point. 


is much smaller than the h x 15 kHz gap to higher bands. 
Consequently, there exists a regime of forces where transi¬ 
tions to higher bands are suppressed and the system is well- 
approximated by a two-band model [19]. 

We probe the lattice geometry with a nearly pure Bose- 
Einstein condensate of ^^Rb, which is initially loaded into the 
lowest band at quasimomentum q = F, the center of the Bril- 
louin zone (BZ) (Fig. 2B). To move the atoms in reciprocal 
space, we linearly sweep the frequency of the beams to uni- 
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FIG. 3. Measuring mixing angles Oq at different quasimomenta. 

(A) Because of the three-fold-rotational symmetry of the honeycomb 
lattice, a path from F to r+3G is equivalent to a triangle-shaped path 
with each leg of length |G|, beginning and ending at F. Coloured 
dots correspond to coloured quasimomentum labels in (B). (B) The 
population remaining in the first band after transport to final quasi¬ 
momentum q. Theory lines are a single-particle solution to the dy¬ 
namics using a full lattice potential and including the first six bands 
(solid) and a two-band, tight-binding model (dashed) [19]. The inset 
Bloch sphere depicts the transported state at F (red), F + G (blue), 
and F + 2G (green) in the basis of the cell-periodic Bloch functions 
at F. Error bars represent the standard error of the mean from averag¬ 
ing 9-11 shots, with the exception of q = M+G and q = M+0.9G, 
which show the average of 20 shots. (C) Transport of a Bloch state 
by one reciprocal lattice vector corresponds to a 27r phase shift in 
the real-space wavefunctions of each sublattice site. Projecting the 
combined lattice and gradient potential V(x) along the path shown 
in red onto the v-axis, which is the direction of the applied force in 
the measurements of Fig. 2B and Fig. 3B, highlights the effect of the 
real-space embedding of the honeycomb lattice. Since the distance 
between A (solid circles) and B sites (open circles) is 1/3 the distance 
between sites of the same type, there is a phase difference of 27r/3 
between the real-space wavefunctions of A and B sites, which gives 
rise to the band mixing. 


formly accelerate the lattice, thereby generating a constant in¬ 
ertial force in the lattice frame. By independently controlling 
the frequency sweep rate of two beams (see Fig. 2A), we can 
tune the magnitude and direction of the force and move the 
atoms along arbitrary paths in reciprocal space. 


To verify that we can access the Wilson line regime, where 
the dynamics are governed entirely by geometric effects, we 
transport the atoms from F to different final quasimomenta 
using a variable force |F| and perform band mapping [17] to 
measure the population remaining in the lowest band (Fig. 
2C). For vanishing forces, we recover the adiabatic limit, 
where the population remains in the lowest band. For increas¬ 
ing forces (i and ii in Fig. IB), where the gradient \F\d over 
the distance between A and B sites d is less than the combined 
width 5, the population continuously decreases. However, at 
strong forces (iii in Fig. IB), where \F\d > s, the population 
saturates at a finite value. For example, after transport by one 
reciprocal lattice vector (blue data in Fig. 2C), about one quar¬ 
ter of the atoms remain in the first band, in stark contrast to 
typical Landau-Zener dynamics, where the population van¬ 
ishes for strong forces [22] . 


Theoretically, the population in the first band after the 
strong-force transport directly measures the Wilson line ele- 


ment|FFpl^q| = (T>q|Wr^q|^f) in the basis of the band 
eigenstates. Based on Eq. 2, the saturation value |FFpl^q|^ = 


I (^ql^f) r population after transport to q is a measure 
of the overlap between the cell-periodic Bloch functions of the 
first band |tiq) at F and q. Notably, for the case of transport 
by one reciprocal lattice vector G, the cell-periodic parts \u'^) 
are not identical, despite the unity overlap of the Bloch states 
l^q) at F and F + G. In contrast to the typical Landau-Zener 
case, they are also not orthogonal—hence the finite saturation 
value. 


To corroborate that our experiment measures the Wilson 
line, we transport atoms initially in the ground state at F by up 
to three reciprocal lattice vectors (Fig. 3). The three-fold ro¬ 
tational symmetry of the lattice, combined with the symmetry 
of its v-orbitals, makes the path from F to F-1-3G equivalent to 
the triangular path shown in Fig. 3A, such that the overlap be¬ 
tween cell-periodic components of the Bloch wavefunctions 
at the two endpoints is unity (see Eq. 2). Correspondingly, we 
expect to recover all population in the lowest band after trans¬ 
port from F to F + 3G. This prediction is confirmed in Fig. 
3B, where we plot the population remaining in the first band 
after transport to final quasimomentum q. The data are well 
described by a tight-binding model that takes into account the 
relative phase between orbitals on A and B sites of the lattice 
due to the Wilson line Wr^q = Physically, this can be 
understood by assuming that the real-space wavefunction sim¬ 
ply accumulates a position-dependent phase when the strong 
force F = /i(q —Q)/t is applied for a short time t (Fig. 3C). 
Notably, the result depends crucially on the real-space embed¬ 
ding of the lattice and would be different in, e.g., a brick-wall 
incarnation [23] of the same tight-binding model. Discrep¬ 
ancies from the tight-binding model result from population 
transfer to higher bands [19]. 

As the Wilson line enables a comparison of the cell- 
periodic Bloch functions at any two quasimomenta (Eq. 2), 
it can in principle be applied to fully reconstruct these states 
throughout reciprocal space. To this end, it is convenient to 
represent the state |tXq) at quasimomentum q in the basis of 
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cell-periodic Bloch functions |1) = |i4q) and |2) = \uq) at a 
fixed reference quasimomentum Q as 

lO=cos^|l) + sin^e*^o|2). (3) 

Mapping out the geometric structure of the lowest band there¬ 
fore amounts to obtaining 6>q and 0q, which parametrize the 
amplitude and phase of the superposition between the ref¬ 
erence Bloch states, for each quasimomentum q [24, 25]. 
Whereas the total phase of |i4q) is gauge dependent, i.e., it 
can be chosen for each q, the relative phase 0q is fixed for all 
q once the basis states |1) and |2) are fixed. Throughout this 
work, we choose the basis states at reference point Q = T. 

In this framework, the population measurements in Fig. 
3B constitute a reconstruction of the mixing angle 0^ = 
2 arccos | fFpl^q |. This can be visualized as the rotation of 
a pseudospin on a Bloch sphere, where the north (south) pole 
represents |1) (|2)). As a function of quasimomentum q, the 
angle 6>q winds by 27r/3 per reciprocal lattice vector (see inset 
of Fig. 3B). 

To obtain the relative phase 0q, which is directly connected 
to the Wilson line via 0q = Arg[IFQ^_^q] - Arg[I^Q^_^q] [19], 
we perform a procedure analogous to Ramsey or Stiickelberg 
interferometry [26, 27]. We initialize atoms in the lowest band 
at F — G and rapidly transport them by one reciprocal lattice 
vector to prepare a superposition of band eigenstates at the ref¬ 
erence point F (i in Fig. 4A). We then hold the atoms at F for 
a variable time (ii), during which the phase of the superposi¬ 
tion state precesses at a frequency set by the energy difference 
between the bands at F. Following this preparation sequence, 
we rapidly transport the superposition state to a final quasi¬ 
momentum q^, lying at angular coordinate a on a circle of 
radius |G| centered at F. Measuring the population of the first 
band as a function of hold time yields an interference fringe 
that reveals the relative phase [19]. 

We observe quantized jumps of tt in the phase of the inter¬ 
ference fringe each time a is swept through a Dirac point, i.e., 
every 60° (blue circles in Fig. 4B) [28, 29]. The binary nature 
of the phases is a consequence of the degeneracy between A 
and B sites, which dictates that the band eigenstates at each 
quasimomentum be an equal superposition of states |4>q) and 
|4>q) on the A and B sublattices [19]. Therefore, on the Bloch 
sphere, the pseudospin is constrained to rotate on a meridian 
about an axis whose poles represent the corresponding cell- 
periodic functions kq) and ) (inset of Fig. 3B). When 
we remove this constraint by introducing an energy offset be¬ 
tween A and B sites [19, 30], we observe smoothly varying 
phases that are always less than tt (red circles in Fig. 4B). 
The dependence of the phase on angle a indicates both the 
breaking of inversion symmetry and the preservation of the 
three-fold rotational symmetry of the lattice. 

Apart from reconstructing the cell-periodic Bloch func¬ 
tions, our method also provides access to eigenvalues of 
Wilson-Zak loops, Wq^q+c. which is essential for deter¬ 
mining various topological invariants [7-9]. To this end, we 
split the Wilson-Zak-loop matrix into a global phase factor, 
which can be measured by extending previous methods [13- 
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FIG. 4. Measuring relative phases 0q at different quasimomenta. 

(A) Schematic of the interferometric sequence in the extended BZ 
scheme (left) and the corresponding rotation on the Bloch sphere 
(right). To create a superposition state, atoms initially in the low¬ 
est eigenstate at T — G are rapidly transported to T (i). The phase 
of the superposition state is controlled by varying the hold time at T 
(ii). After the state preparation, the atoms are transported to a final 
quasimomentum q^, which is parametrized by the angle a and lies 
on a circle of radius |G| centered at T (iii). (B) Phases referenced 
to a=180° for the lattice with AB-site degeneracy (blue) and AB-site 
offset (red). Data in blue have been offset by +120° for visual clarity. 
Dashed lines are a two-band, tight-binding calculation with A/ J=0 
(blue) and A/J = 3.1 (red), where J — hx 500(10) Hz. Error bars 
indicate fit errors. 


15, 31], and an SU{2) matrix with eigenvalues Us¬ 

ing the data from Fig. 3B and Fig. 4B, we reconstruct the 
eigenvalues for a loop transporting from F to F + G, up to 
multiples of tt [19]. We find the eigenvalue phases to be 
^ = 1.03(2)7r/3, in good agreement with the value of 7r/3 
predicted from the two-band model. Remarkably, we mea¬ 
sure the same eigenvalues even when the band eigenstates are 
modified by an energy offset between A and B sites [19]. This 
invariance is a direct consequence of the real-space represen¬ 
tation of the Wilson-Zak loop, Wp^r+G = (see Eq. 2 
and [19]). Because the Wilson-Zak loop depends only on the 
position operator f, the eigenvalues are determined solely by 
the physical locations of the lattice sites, which are unchanged 
by the energy offset. 

Our versatile approach to measuring Wilson lines only 
employs standard techniques that are broadly applicable 
in ultracold atom experiments, and can be extended to 
higher numbers of bands by adopting ideas from quantum 
process tomography [32]. Provided that the relevant local 
Hilbert space is identical for all quasimomenta, our method 
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provides a complete map of the eigenstates over the BZ, 
giving access to the Berry curvature and Chem number. 
When this is not the case, the same techniques enable 
the reconstruction of Wilson-Zak loop eigenvalues, which 
directly probe the geometry of the Wannier functions [9] and, 
therefore, the polarization of the system [21, 33]. Conse¬ 
quently, these eigenvalues can reveal the topology of bands 
with path-dependent and non-Abelian Wilson lines [9, 11]. 
Such systems can be realized in cold atom experiments 
by periodically modulating the lattice [12, 30, 34-36] to 
create a quasimomentum-dependent admixture of additional 
bands [37] or coupling between internal states [38-41]. 
Moreover, the addition of spin-orbit coupling [42] would 
enable the investigation of the Z 2 invariant characterizing 
time-reversal invariant topological insulators [34, 41, 43-46]. 
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SUPPORTING ONLINE MATERIAL 


SI. THEORETICAL BACKGROUND 


In these supplements, we provide theoretical background 
on our measurement of Wilson lines in a honeycomb lattice 
and additional experimental details. In SI A, we derive the 
equation of motion for a single particle in a lattice in the pres¬ 
ence of a constant force, relate it to the Wilson line, and apply 
our results to the honeycomb lattice. We then describe our 
numerical calculations of the Wilson line in SIB. In SIC, we 
relate the Wilson line to the cell-periodic Bloch functions of 
the lattice and the position operator. We then check the valid¬ 
ity of assuming completeness of the bands in the experiment 
in SID. Next, we discuss general gauge freedom in quantum 
mechanics and its application to our measurements in SIE. 
We then present a decomposition of the f/(2) Wilson line into 
a U (1) part and an SU ( 2 ) part in SIF. We conclude the the¬ 
oretical section of the supplements by describing the recon¬ 
struction procedure for the S'f/(2) eigenvalues of the Wilson 
line in SIG. Experimental methods, including data analysis 
techniques, are given in SII. 
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A. Dynamics in the combined lattice and gradient potential 


Here, we derive the equations of motion for a particle in a 
lattice in the presence of a constant force F. The Hamiltonian 
of the lattice can be written as 

F = (S.l) 

q,n 

where is the energy of the nth band at quasimomentum 
q and |^q) are the Bloch states. The Bloch states can be ex¬ 
pressed as l^q) = where jn^) are the cell-periodic 

Bloch functions and f is the position operator. 


Adding a constant force F to the system results in the 
Schroedinger equation: 

idt\m) = iH-F-r)m)), (S.2) 

where we have taken h = 1. We assume the initial state is 
localized in reciprocal space at quasimomentum qo such that 

|V’(0))=5]a"(0)|$"„) (S.3) 


where |a’^( 0 )|^ gives the population in the nth band at time 
t = 0. Making the ansatz 

= (S.4) 

n 

q(0 = qo + Fi (S.5) 

leads to the following equations of motion for a two-band sys¬ 
tem: 




^a(t 


2,1 

q(^) 


q(^) 


^q(t) 

^^(0 - ^ 



faHt)\ 

Wit)) 

(S.6) 


where 


f-n,n' _ A n,n' j-y 

^q(t) - ^q(t) • ^ 

= (S-7) 

and 

AqS'=*K|Vql<)lq=q(t)- (S-S) 

The quantity defines an intra-band (n = n') and an 

inter-band (n 7 ^ n') Berry connection. From Eq. S. 6 , we 
see that the inter-band Berry connection drives transitions be¬ 
tween the different bands. 
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i. The limit of infinite force 


In the limit of an infinite force, the energy terms on the di¬ 
agonal are negligible compared to the geometric terms . 
In this case, Eq. S.6 reduces to 


idt 


Wit)] 




C(J 


Wit)] 


(S.9) 


Defining as the matrix with elements , the evolution 
is given by 


\W))=Texp[i j di|q(t)]|V»(0)) 

= W|t/)(0)) (S.IO) 


Using Eq. S.5 to change the variable of integration from time 
to quasimomentum space recovers Eq. 2 of the main text: 


Wq(o)^q(t) = Pexp[i^dqAq], (S.ll) 

where the path-ordered {V) integral runs over the path C in re¬ 
ciprocal space from q(0) to q(t) and Aq is the Wilczek-Zee 
matrix with elements defined in Eq. S.8. Therefore, the evo¬ 
lution of the system is described by an operator W, which is 
the Wilson line [SI]. Path-ordering {V) is necessary because 
the matrices generally do not commute for all quasimomenta 
along the path. 


a. Dynamics in the tight-binding honeycomb lattice 



FIG. SI. The honeycomb lattice in real space and re¬ 
ciprocal space. A, The real-space honeycomb lattice com¬ 
prises triangular sublattices A (solid circles) and B (open cir¬ 
cles) with nearest-neighbour hopping vectors di. B, Pseudospin 
S(q) = (sin cos 0q, sin sin (/)q, cos ^q) for a lattice with AB- 
site degeneracy in a basis formed by the cell-periodic Bloch states at 
reference quasimomentum Q = T, as labelled. The components Sx 
and Sy are indicated by white arrows, with the length of the arrow 
representing the magnitude of the Sx{y) component; Sz is illustrated 
by the color map, with red (blue) indicating > 0 {Sz < 0). 

We now apply the results from the previous section to the 


specific case of the tight-binding model of the honeycomb lat¬ 
tice. The honeycomb lattice may be decomposed into two 
triangular sublattices composed of A and B sites and coupled 
via nearest-neighbour lattice vectors di with hopping ampli¬ 
tude J (see Pig. SI A). We begin by defining the states |^q) 
and l^q ) of the A and B sites via states |u;rA(B)) localized on 
the A {B) sites as 

= eWnt) (S.12) 

(S.13) 

where N denotes the number of lattice sites. In the tight- 
binding limit, where the localized states are eigen¬ 

states of the position operator, the cell-periodic Bloch func¬ 
tions are independent of quasimomentum and can be 

physically understood as superpositions of the localized wave- 
functions |r^rA(B)) ^11 ^ 

In the basis of |^q) and |^q ), the Hamiltonian describing 
the two lowest bands of the honeycomb lattice is [S2] 

^tb(q) = 

where A is an energy offset between the sublattices and 

iq = \tq\ 

= - (S.15) 

When the A and B sites are degenerate at A = 0, this Hamil¬ 
tonian is diagonalized by the eigenstates 

l^q) = y|(l^q)-e'’'‘'l^q)) (S. 16 ) 

l^q) = ;^(l^q)+e'’'‘’l^q»- (S-17) 

The corresponding eigenenergies are 

= - |tq| (S.18) 

El=\tq\. 


Hi. Elements of the Wilczek-Zee connection for A = 0 


To calculate the connections in the Wilczek-Zee matrix, we 
note that 

l<) = ^i\<) - eWu^)) = (S.19) 

\ul) = E(^\u^)+eWu^))=e-W^l). 


(S.20) 
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The intra-band (Berry) connections are then 

= (S.21) 

and the inter-band connections are 

Af = A2 -i = (S.22) 


B. Numerical calculation of Wilson lines 

Using the equations of motion derived in Eq. S.6, solving 
for the Wilson line transporting from quasimomentum q(t = 
0) to quasimomentum q(t = T) yields 

Wq(0)^q(r) (S.23) 

where H (t) is the matrix on the right-hand side of Eq. S.6. 

To numerically calculate the Wilson line, we compute the 
Trotter product of n time-independent matrices evaluated at 
discrete time-steps of size At = Tln: 

Wq(0)^q(T) (S.24) 

n 

i=i 

In our calculations, we use several hundred time-steps, de¬ 
pending on the length of the path. 

In the tight-binding model, we have analytical expres¬ 
sions for the eigenenergies (Eq. S.18) and the Berry connec¬ 
tions (Eqs. S.21 and S.22); it is therefore straight-forward to 
calculate H{t) for any time t. In the ab-initio calculation, 
obtaining the Berry connections requires a particular gauge- 
choice to ensure that the cell-periodic Bloch functions are nu¬ 
merically differentiable. Similar to the approach in Ref. [S3], 
we choose our gauge such that the Bloch functions of the 5- 
bands are entirely real on a lattice site and the Bloch functions 
of the p-bands are entirely imaginary midway between neigh¬ 
boring lattice sites. This gauge-choice allows us to numeri¬ 
cally differentiate the cell-periodic Bloch functions along the 
Qx direction, which is sufficient for our experiments. 


C. Wilson lines as projectors 

To derive the relation between the Wilson line and the cell- 
periodic Bloch functions, we discretize the Wilson line in 
Eq. S.ll by dividing the path from q(0) to q(t) into N in¬ 
finitesimal segments qi,q 2 ---qN- The Wilson line can then 
be expressed as a sequence of path-ordered products of pro¬ 
jectors ^(q) = |^q)(^q| [S4] with elements 

N 

= ««l n (S.25) 

i=l 


where Af is the number of bands. In the experiment, Af = 2. 

When the states \u^) form a complete basis over the 
Hilbert space the projectors are trivial, i.e. ^(q) = 
l^q)('^ql = 1. In this case, the Wilson line elements 
reduce to the overlap of the cell-periodic Bloch functions: 

W^qT->q. = KJO- (S- 26 ) 

Therefore, the Wilson line provides a way of comparing the 

cell-periodic Bloch states |r^q) at any two points in momentum 
space. In terms of the Bloch states |T>q) = the 

Wilson line elements can equivalently be expressed as 

^(^n^|g,Aq.r|^™) (s_ 27 ) 

where Aq = — qi is the change in quasimomentum. 

Therefore, expressed in terms of the real-space position op¬ 
erator f, the Wilson line transporting a state from an initial 
quasimomentum q^ to final quasimomentum q/ by Aq is 

(S.28) 


D. Validity of the two-band tight-binding approximation 

In this section, we ascertain the validity of applying 
Eq. S.26 to our system, which only holds when the two-lowest 
band eigenstates span the same Hilbert space at all quasi¬ 
momenta. To this end, we use the numerically-calculated 
band eigenstates of the full optical lattice potential to calcu¬ 
late the elements of the Wilczek-Zee connection for 

= 1,2. We compare the corresponding Wilson lines 
at various lattice depths to the Wilson line obtained from a 
tight-binding calculation. 

We plot the population in the first band, | | ^ at dif¬ 

ferent lattice depths for a path from T to T -f 3G in Eig. S2. 
While strong deviations occur for shallow lattices, only minor 
differences are visible at the lattice depth of b.2Er used in the 
experiment. Therefore, the experiment should, in principle, 
be well-described by the tight-binding formalism. 

However, in Eig. 3 of the main text, the population in the 
lowest band returns to only 90% and not unity at T -|- 3G. 
This is a result of transfer into higher bands. Due to the pres¬ 
ence of higher bands, we can not exactly realize the two-band 
Wilson lines plotted in Eig. S2, which would require reach¬ 
ing the infinite gradient limit for the two lowest bands while 
remaining adiabatic with respect to higher bands. Experi¬ 
mentally, the choice of gradient strength is a compromise be¬ 
tween realizing dynamics that are fast compared to the energy 
scale of the lowest two bands and minimizing excitations into 
higher bands. We could, in principle, realize the two-band 
Wilson line regime more precisely by increasing the lattice 
depth, which decreases the combined width 5 of the lowest 
two bands and increases the energy scale between the lowest 
two bands and higher bands. However, in the current work, 
the lattice depth was limited by the bandmapping technique. 
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As the lattice depth is increased, it becomes more difficult for 
the bandmapping process to remain adiabatic with respect to 

5 . 



Final quasimomentum 

FIG. S2. Validity of the two-band tight-binding approximation. 

Lowest band population, | | ^ for a two-band Wilson line using 

a full band structure calculation for lattice depths S.lEr, 3Er, 2Er, 
and \Er (decreasing lightness of blue). The Wilson line at b.2Er 
is nearly identical to the Wilson line obtained from a tight-binding 
calculation (orange). 


E. Gauge Freedom in Wilson lines 

Due to its linear structure, quantum mechanics contains an 
inherent gauge freedom: if 1 7 /^) is an eigenstate of an operator, 
then so is (j) ^ M. Consequently, even non-degenerate 

eigenstates are defined uniquely only up to such a phase fac¬ 
tor. If the Hamiltonian depends on a parameter R as ^(R), 
then the corresponding phase 0(R) can be chosen indepen¬ 
dently for every value of R. This is called a local (7(1) gauge 
freedom, in reminiscence of the situation in quantum electro¬ 
dynamics or gauge theories [S5]. 

Due to this gauge freedom, the Berry phase for adiabatic 
evolution in a non-degenerate band is only well-defined for 
closed loops, since, for non-closed paths, the phase factor in 
front of the final state depends on the choice of basis, i.e., the 
gauge choice, at the final point. In the situation of two, every¬ 
where non-degenerate eigenstates 11) and |2), the gauge group 
enlarges to (7(1) x (7(1), since the phase of each eigenstate 
can be chosen independently. However, if the two eigenstates 
are degenerate with eigenenergy E, then any normalized su¬ 
perposition |T^) = a\l) -h 13\2) is also an eigenstate with the 
same eigenvalue: 


^|^) = ^(a|l)+^|2)) 

= aH\l) ^ f3H\2) 

= E\^) (S.29) 

The freedom in choosing the basis states is therefore now en¬ 
larged. For each normalized superposition state |T^), there ex¬ 
ists a corresponding orthogonal superposition state such 
that |T^) and form an orthonormal basis. Consequently, 


the gauge freedom is U(2). 

The Hamiltonian is, however, not the only observable. If 
there exists another observable that is conserved by all ^(R) 
and can distinguish between the states |1) and |2), such as 
spin or parity, then this observable defines a new, more con¬ 
strained basis. Assuming the basis states to be eigenstates of 
both the Hamiltonian and the observable, the gauge freedom 
is reduced again to (7(1) x (7(1), even for degenerate eigenen- 
ergies. This is precisely the case in the strong-force limit of 
the experiment: even though the bands appear degenerate dur¬ 
ing the evolution, i.e., the difference in eigenenergies is neg¬ 
ligible during the dynamics, the bandmapping procedure can 
nonetheless distinguish between the two bands. In such a ba¬ 
sis, the absolute values of the Wilson line elements are well- 
defined and can be observed even for open lines. 


/. The pseudospin representation and the Wilson line elements 


We now explicitly discuss the gauge-invariance of the terms 
6>q and used in the pseudospin representation of the cell- 
periodic Bloch state. The components of the pseudospin ref¬ 
erenced to Q = r is plotted in Fig. SIB for a lattice with 
AB-site degeneracy. 

The polar angle 6>q is obtained from the quantity 





^ 6>q = 2arccos |(rtq|l)| 
= 2arccos | fCg^-^q 


(S.30) 


where we have used Eq. S.26 to relate the overlap of the cell- 
periodic Bloch functions to the Wilson line elements in the last 
line. The (7(1) x (7(1) gauge-freedom on the phase of |?Xq) 
or the reference states |1) and |2) does not affect the absolute 
value of the overlap. Hence, | \ and, consequently, 6>q 

is a gauge-invariant quantity. 

The relative phase 0q can be expressed as 

(/)q = Arg[(ui |1)] - Arg[(ui |2)] 

= - Arg[W^\^] (S.31) 

That is, we access the relative phase between the basis states 
by obtaining the difference of phases between the Wilson line 
elements. While the gauge-freedom of |rtq) is cancelled out 
by taking the difference of the argument of Wilson line ele¬ 
ments, the gauge-freedom on the reference states remains. A 
different choice of the phase for |1) and |2) changes the value 
of 0q. However, to obtain a gauge-invariant quantity, we can 
compare (/)q and 0q/ at quasimomenta q and q'. Taking the 
difference between the relative phases at the two quasimo¬ 
menta cancels out the gauge-freedom of the reference states 
if the same reference states are used to define the cell-periodic 
Bloch function at both quasimomenta. Explicitly, the gauge- 
invariant quantity measured in the experiment is 


0q : 


(S.32) 
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R Decomposition of the Wilson line into U(l) and SU(2) 


Here, we decompose the U(2) Wilson line into an SU(2) 
matrix multiplied by a global U(l) phase. We furthermore 
relate the U(l) phase to the sum of the Berry phase of the first 
and second band, noting, however, that the U(l) part might 
not show up in certain bases [S6], and that our experiment 
measures only the SU(2) part of the Wilson line. 

We begin by writing the Wilczek-Zee connection matrix 
Aq in a general form as 


A 


q ~ 




(S.33) 


To simplify notation, we henceforth suppress the q subscript 
and note that all elements are still to be understood as being 
q-dependent. Next, we decompose the matrix as 


A = 




0 



where = 1. The eigenvalues of this 

matrix depend on the absolute values of the Wilson line terms 
and the phase of and are given by 

± + ini[W^ii]2 (S.38) 

The absolute values are directly measured via the popula¬ 
tion remaining in the first band after transport of the lowest 
band eigenstate. Although the interferometric sequence re¬ 
veals only the difference between the phases of elements 
and we can extract the phase of by invoking uni- 
tarity and the ’’backtracking” condition of Wilson lines. It can 
be shown that, for generic quasimomenta q and Q, Wg^q = 
Wq^Q [S9]. Consequently, Wq^QWg^q = 1, such that 
going forward and back along the same path results in no 
transformation of the state vector. 

In the case of the Wilson-Zak loops in the experiment, the 
relevant relation is 


Wr^r+G = ^r^r-G (S-39) 


•— + Asu{2) (S.34) 

Noting that A^^i) is proportional to the identity matrix and 
therefore commutes with Asui 2 ) and with itself for all quasi¬ 
momenta, the Wilson line Wg^q transporting a state from 
initial quasimomentum Q to final quasimomentum q can be 
expressed as: 

Wq^q = 

= Ve^ ^qA«7(i)+As«7(2) 


where we have used that Wr^r-G = Wr+G^r which, as¬ 
suming a periodic gauge choice, follows from the periodicity 
of the BZ. If the phase (j) of the oscillation after transport from 
r to qo, = r H- G is given by 

^ = + Arg[VF/^qJ - Arg[WP_^J (S.40) 

where (j)p is a gauge-dependent quantity that results from the 
initial transport to prepare the superposition state, then the 
phase 0' of the oscillation after transport from T to T — G 
is 


_ gi dqA[7(i) dqAsu(2) (S.35) 


(j)' _ (j)^ — Arg[kFpl^p^Q] — Arg[kFp^p^Q] + tt (S.41) 


where C denotes the path taken from Q to q. This gives the de¬ 
composition of the (2) Wilson line into a f/(l) global phase 
multiplied by a path-ordered SU(2) matrix. Furthermore, the 
global U{1) phase is given by 


Therefore, taking the difference between the two oscillation 
phases extracts Arg[IRpl^p^Q] as 

- 0' - TT = 2Arg[W]i^r+G] (S.42) 


Jc ^ 2 


(S.36) 


where is the adiabatic phase acquired in the first band and 
02 is the adiabatic phase acquired in the second band. For 
closed loops, 01 and 02 are the Berry phases, which can be 
used to formulate the Chern number of the system [S7] . For 
paths closed only by a reciprocal lattice vector, 0i and 02 are 
instead the Zak phases [S8]. 


G. Reconstruction of SU(2) Wilson-Zak loops 


Generically, an SU(2) matrix can be expressed as 


/ IRii IR12\ 


(S.37) 


Note that there is an ambiguity in choosing ±7r when relating 
0 ' to 0 . This results in a global U{1) phase shift of tt in the 
eigenvalue phases. However, the phase difference between the 
eigenvalues is unaffected, which is sufficient to reconstruct, 
e.g., the Z 2 invariant. 

When the A and B sites of the lattice are degenerate, 
applying Eq.S.42 to the phase of oscillations for a = 0 
and a = 180 in the interferometric sequence (Fig. 4B 
of main text) yields Arg[H7pl^p^Q] =0.03(7) rad. Com¬ 
bined with the direct transport data from F to F - 1 - G 
(Fig. 3B of main text), which gives |H^rl^p^Q | =0.47(2) and 

|kkr^p+G I = — IkFpl^p+Gl^ =0.88(1), we obtain eigen¬ 

values exp[±il.03(2)7r/3]. 

Similarly, in the lattice with AB-site offset, data for 
a = 0 and a = 180 from Fig. 4B in the main text yields 
Arg[H7pl^P^Q] =-0.76(6) rad. We measure the absolute values 
by transporting atoms initialized at F in the lowest eigenstate 
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FIG. S3. Reaching the Wilson line regime in a lattice with AB-site 
offset. The population remaining in the first band after transport at 
different forces from T to T + G. The data agrees reasonably well 
with a two-level, tight-binding theory (dashed line) that approaches 
the Wilson line regime (thick shaded line) at large forces. We at¬ 
tribute the discrepancy to the two-level model at larger forces to 
transfer to higher bands. To calculate the SU(2) eigenvalues, we use 
the population at |F| d/e = 5. The inset depicts the transport path. 
Error bars represent the standard error of the mean from ten shots per 
data point. 


to r + G with increasing force. The remaining population 
in the first band is shown in Fig. S3. The eventual saturation 
of population transfer indicates the geometric nature of the 
transfer. At |F| dje — 5, we obtain | =0.63(3) and 

|VFp^p^q| = — |VFrl^p^Q|^=0.77(2). The eigenvalues 

of this Wilson-Zak loop are then exp[±il.04(4)7r/3]. 


SII. EXPERIMENTAL METHODS 


A. The optical potential of the honeycomb lattice 


The total potential resulting from interfering three beams of 
variable polarization at 120° angles can be decomposed into 
the sum of its out-of-plane (s-) and in-plane (p-) components 
as 


V{x,y) =V^{x,y) + VP{x,y) 



i=l 



(S.43) 


where is the ac Stark shift produced by the 

s(p)-component, is the wave-vector with wavenumber 
= \ki\^ and ai is the phase between s- and p-polarization 
components of beam i. 

In the lattice with AB-site degeneracy, all three beams have 
equal intensity and are purely 5-polarized. In this case, the 


expression for the total potential reduces to 
V{x,y) = Vb(^2cos(y3/cLx) 

+ 4cos(2^^)cos(5^)+ 3 ) (S.44) 

where Vo = = Vi. 

The resulting honeycomb potential contains two non¬ 
equivalent lattice sites (A,B) per unit cell, as shown in Fig. 
S4. Consequently, the two lowest bands, which correspond to 
the v-orbitals on the A and B sites, touch at Dirac points and 
are strongly coupled to each other. Coupling to the next higher 
bands, however, can mostly be neglected due to the large en¬ 
ergy gap to the p-orbitals. Hence, our experimental system is 
well-approximated by a two-band model. 


i. Breaking AB-site degeneracy 

To introduce an energy-offset between the A and B sites 
and maintain isotropic tunnelling, we constrain the polariza¬ 
tions of each beam to have the same composition of s- and 
p-polarizations. The potentials arising from the interference 
of the s- and p- components of the three beams, which are 
shown separately in Fig. S4, have the form 


V^{x,y) = Vbcos^6>[ 2 cos{V3kLx) 


VSkLX MlV. , „ 
+ 4 cos(—^—) cos(^—) + 3 


and 


y) = — Vo sin^ 0 ( cos( 


s/SkLX 3kLy 


+ 


(S.45) 


■ <^ 32 ) 


V^kLX 3kLy , , . 

+ COs(- - - - -h 0^13 + <^32) 

+ cos{V3kLX + ais) — 3 ) (S.46) 


where 0^32 = 0^3 — 0^2 and 0^13 = o^i — 0^3 and 0 parametrizes 
the composition of 5- and p-polarizations, i.e., for 0=0, the 
light is purely s-polarized and for = 7r/2, the light is purely 
p-polarized. Furthermore, in defining the same Vq for the s- 
and p-polarizations, we have neglected the state-dependence 
of the dipole potential, which is valid in our case of far- 
detuned light. 

By choosing the phase of each beam, we can shift the p- 
polarized potential relative to the s-polarized potential. When 
<^32 = <^13 = 0 , the minima of the p-polarized potential and 
the maximima of the s-polarized potential coincide with the 
A and B sites (Fig. S4A). Subsequently, atoms experience the 
same ac Stark shift in either an A or B site. However, by set¬ 
ting 0^32 = 27r/3 and ais = 27r/3, the s-polarized potential is 
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FIG. S4. The honeycomb lattice potential, a, A lattice with degenerate A and B sites formed by beams of polarization angle ^ = tt/G and 
0^32 = tti 3 = 0. Top: 2D plots of the s-polarized (left) and p-polarized (right) components of the potential. In the s(p)-polarized potential, 
the A (gray circles) and B sites (green circles) are both located at the potential minima (maxima). Therefore, there is no energy offset in the 
total potential, which is the sum of the two polarization components at the A and B sites. Bottom: A cross-cut of the potential through the 
dashed line in the 2D plots. The s(p)-potential is in blue (red) and the total potential is in green, b, Same as in (a), but for 0^32 = 27r/3 and 
ai 3 = 27r/3. With the appropriate phase shift between the polarization components of the beams, the A sites are located at the minima while 
the B sites are located at the maxima of the p-polarized component of the potential. Consequently, there is an energy offset between A and B 
sites in the total potential. 


shifted such that the potential maxima occur on A sites while 
the potential minima occur on B-sites (Fig. S4B). 


a. Implementation of the honeycomb lattice 

For the lattice with AB-site degeneracy, the polarization of 
the three beams is set by polarizing beam splitters. We have 
verified in previous work [SIO] that this results in sufficiently 
pure s-polarizations. 

To introduce an AB-site offset, we tune the polarizations of 
the beams by using a half- and a quarter-waveplate in the paths 
of two beams and only a half-waveplate in the path of the 
third beam, which does not require a phase shift between its 
s- and p-polarized components. After setting the waveplates, 
we ensure that the polarization composition of each beam is 
approximately equal by taking time-of-fiight (TOF) images of 
the BEC after sudden release from the lattice. An unequal 
polarization composition between the beams results in an im¬ 
balance in the Bragg peaks. We then check the dispersion re¬ 
lation through the Ramsey-like interferometric procedure de¬ 
scribed in Sec. SI G. To quantitatively assess the amount of 
AB-site offset, we fit the measured dispersion relation to a 
tight-binding model (see Fig. S7). 


B. Preparation scheme 

The evaporative cooling of ^^Rb atoms in the 
|F = I, mi? = I) state to quantum degeneracy is initi¬ 
ated in a plugged quadrupole trap and completed in a 
crossed-beam dipole trap of wavelength 1064 nm. At the 
end of the cooling process, we have approximately 4x10^ 
atoms in the BEC. The atoms are adiabatically loaded 
into a honeycomb optical lattice of depth 5.2(I)E^^ in 100 
ms. During the experimental sequence, the combined trap 
frequencies of the blue-detuned lattice and dipole potential 
are ujz = 118(9) Hz and ujxy = 16(1) Hz. We obtain these 
frequencies by measuring the oscillation frequency of the 
center-of-mass motion of the BEC after a perturbation of the 
trapping potential. We neglect the effect of the dipole trap 
since the dynamics of our experiment (on the order of 500 /is) 
is much shorter than the inverse dipole trap frequencies. 


C. Lattice acceleration 

To transport the atoms in reciprocal space, we generate a 
constant inertial force in the lattice frame by uniformly ac¬ 
celerating the lattice. An acceleration of in 

the propagation direction of beam i is accomplished via a 
linear sweep of the frequency shift Ui = Sui/(27r) at a rate 
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^. Individual control over the frequency sweep rate of two 
beams enables lattice acceleration of variable magnitude and 
direction. Thus, we can move the atoms along arbitrary paths 
in reciprocal space. 



FIG. S5. Raw data of band mapped atoms at F. a, The extended 
zone scheme showing the first (hexagon) and second (triangles) BZs 
is overlaid on a raw image of the band mapped atoms at F. High- 
symmetry points F, at the center of the first BZ, and M, at the edge 
of the first BZ, are labelled. Non-equivalent Dirac points K (K') 
are depicted by solid (open) orange circles at the corners of the first 
BZ. b, Analysis ROIs. The atom number in the first (second) zone is 
obtained by summing the pixel values within the red (blue) circle(s). 
We additionally take the mean of the pixel values in the yellow ring 
and, with the exception of the interferometric data, subtract this value 
as background from the pixel sum of the first zone atoms. 


band at F are unstable due to interaction effects and move 
to other quasimomenta. An additional background subtrac¬ 
tion would therefore underestimate atoms in the second band 
by counting atoms that have decayed from F as background. 
This analysis method was used to obtain the population in the 
lowest band shown in Figs. 2 and 3 in the main text. 

In contrast, for the oscillation data in Fig. 4B of the main 
text, we do not perform a background subtraction forni. Due 
to the long hold times, the cloud heats and disperses in re¬ 
ciprocal space, increasing the background value. Therefore, 
background subtraction would lead to an underestimation of 
atoms in the first zone. However, an overall offset in the first 
zone population does not affect the phase of the oscillation, 
which is the relevant quantity. 

To check systematic errors due to our selection of ROIs, 
we analyse a single dataset of population transfer vs. force 
magnitude after transport by one reciprocal lattice vector (blue 
data in Fig. 2B of main text) using different ROIs. We eval¬ 
uate both the effect of the ROI size using a fixed background 
subtraction ring and the effect of the background subtraction 
ring using a fixed ROI size. Using the same ROI size for first 
and second band atoms and restricting the ROI size such that 
it does not overlap with the background subtraction ring yields 
consistent results with deviations on the order of ±5%. 

We use the same ROIs and, when applicable, background 
subtraction ring for all datasets. 


D. Detection 

The detection procedure begins with a linear ramp-down of 
the lattice intensity in 800 fis to band map the atoms. We then 
use absorption imaging to detect the atoms after 9 ms TOF. 
Due to the short TOF, the resulting image is a convolution of 
the insitu cloud size and the quasimomentum. Nonetheless, 
for quasimomenta near the center of the first BZ, we can eas¬ 
ily distinguish between atoms in the first and second BZs, as 
shown in Fig. S5A. In contrast, it is difficult to differentiate 
between first and second zone atoms at the edges of the BZ. 
For these quasimomenta, we add an additional adiabatic seg¬ 
ment to the sequence to push the atoms away from the edge, 
toward the center of the BZ, before bandmapping. 


i. Data Analysis 

We sum the pixel values of the atoms in the first zone ni and 
the pixel values of the atoms in the second zone n 2 to obtain 
the fraction of atoms in the lowest band, ni/{ni + 712 ). Since 
we wish to count atoms localized at specific quasimomenta, 
we specify regions of interest (ROIs), which are depicted in 
Fig.SSB for atoms at F. We sum pixel values within the red 
circle to obtain ni and sum the pixel values within the six blue 
circles to obtain 77 - 2 . 

For a quantitatively accurate fraction in the lowest band, we 
subtract the mean pixel value of the shaded yellow region from 
rii to account for the hot background atoms. We do not sub¬ 
tract an additional background for 722 since atoms in the upper 


E. Fitting the interference fringe 

The population in the first band Pi{t) resulting from the 
Ramsey-like interferometric sequence described in the main 
text oscillates as a function of hold time and is given by: 

Pi (t) = Co + Aocos{et + (/)) (S.47) 

where Co is a constant offset, Aq parametrizes the amplitude 
of the oscillation, and the phase 0 is given by Eq. S.40. 

To extract the phase of the interference fringe, we fit the 
population of the lowest band at quasimomentum to an 
empirically chosen function of the form: 

(cos(27r/i + </.)+ yi) + yo (S.48) 

where Aq is the amplitude of the function, to parametrizes 
the decay of the fringe, / gives the frequency, which is de¬ 
termined by the dispersion at the reference quasimomentum, 
and (j) is the phase. The offsets yi and yo interpolate between 
an oscillation with damping of both maximum and minimum 
values (Fig. S6A) and an oscillation with damping of only the 
maximum values (Fig.SbB). 


F. Accessing the dispersion relation 

In addition to probing the band geometry, the interfero¬ 
metric sequence simultaneously reveals the dispersion rela¬ 
tion through the frequency of the oscillation. By varying 






15 


A B 




FIG. S6. Example oscillations from the Ramsey-like interferometric sequence. Both the maximum and minimum values of the interference 
fringe for a = 120° in a lattice with AB-site offset (a) damps with increasing hold time, in contrast to the interference fringe for a = 120° in 
a lattice with AB-site degeneracy (b). Here, only the maximum values damp with increasing hold time. 



the reference quasimomentum Q, we obtain the energy dif¬ 
ference between the lower and upper bands over the entire 
BZ [Sll, S12]. The measured dispersion along the path T-K- 
M-F is shown in Fig. S7. This method is a convenient tool 
for calibrating the lattice depth and quantifying the AB-site 
offset. 


FIG. S7. Mapping the dispersion relation over the BZ. The dis¬ 
persion relation along the high-symmetry paths for lattices with: 
A/ J=0 and depth Vb =0.8F^r (dark blue) or Vb =2.5Er (light blue); 
A/J=3.1 and depth Vb —52Er (red). Theory lines show a full 
band structure calculation (solid) and a best-fit tight-binding model 
(dashed). Error bars indicate fit errors. 


[51] Frank Wilczek and A Zee, “Appearance of gauge structure in 
simple dynamical systems,” Phys. Rev. Lett 52, 2111-2114 
(1984). 

[52] Gordon W. Semenoff, “Condensed-Matter Simulation of a 
Three-Dimensional Anomaly,” Phys. Rev. Lett. 53, 2449- 
2452 (1984). 

[53] W. Kohn, “Analytic Properties of Bloch Waves and Wannier 
Functions,” Phys. Rev. 115, 809-821 (1959). 

[54] Rui Yu, Xiao Liang Qi, Andrei Bernevig, Zhong Fang, and 
Xi Dai, “Equivalent expression of Z 2 topological invariant 
for band insulators using the non-Abelian Berry connection,” 
Phys. Rev. B 84, 075119 (2011). 

[55] M. E. Peskin and D. V. Schroeder, An Introduction to Quantum 
Field Theory (Frontiers in Physics) (Westview Press, 1995). 

[56] A. Zee, “Non-abelian gauge structure in nuclear quadrupole 
resonance,” Phys. Rev. A 38, 1-6 (1988). 


[57] D. J. Thouless, M. Kohmoto, M. P. Nightingale, and M. den 
Nijs, “Quantized Hall Conductance in a Two-Dimensional Pe¬ 
riodic Potential,” Phys. Rev. Lett. 49, 405-408 (1982). 

[58] J. Zak, “Berry’s phase for energy bands in solids,” Phys. Rev. 
Lett. 62, 2747-2750 (1989). 

[59] Yuri Makeenko, Methods of Contemporary Gauge Theory 
(Cambridge University Press, 2005). 

[510] L. Duca, T. Li, M. Reitter, 1. Bloch, M. Schleier-Smith, and 
U. Schneider, “An Aharonov-Bohm interferometer for deter¬ 
mining Bloch band topology,” Science 347, 288-292 (2015). 

[511] A. Zenesini, D. Ciampini, O. Morsch, andE. Arimondo, “Ob¬ 
servation of Stiickelberg oscillations in accelerated optical lat¬ 
tices,” Phys. Rev. A 82, 065601 (2010). 

[512] Sebastian Kling, Tobias Salger, Christopher Grossert, 
and Martin Weitz, “Atomic Bloch-Zener Oscillations and 
Stiickelberg Interferometry in Optical Lattices,” Phys. Rev. 
Lett. 105, 215301 (2010). 








